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FOREWORD 


This  research  was  conducted  by  the  Engineering  Department,  Cambridge 
University,  under  contract  number  ARDSG-27-82,  from  the  European  Research 
Office,  London.  The  studies  were  sponsored  by  the  U.S.  Army  Construction 
Engineering  Research  Laboratory  (USA-CERL)  for  the  Office  of  the  Chief  of 
Engineers  (OCE) .  The  work  was  conducted  under  Project  4A161102AT23,  "Basic 
Research  in  Military  Construction";  Task  A,  "Base/Facility  Development"; 
and  Work  Unit  022,  "Mechanical  Properties  of  Polymers,"  of  which  it  is  one 
of  several  sub-tasks  in  the  research.  The  overall  research  is  designed  to 
provide  a  better  understanding  of  the  micromechanics  of  deformation  and 
rupture  of  polymer  chains  so  that  bulk  material  performance  is  better 
understood  and  acceptable  design  values  for  polymer-based  materials  can  be 
provided  to  engineer  and  architect  designers. 

The  USA-CERL  Principal  Investigator  is  Mr.  Alvin  Smith  of  the  Engineer¬ 
ing  and  Materials  Division  (IM).  Dr.  Robert  Quattrone  is  Chief,  USA-CERL-EM 

COL  Paul  J.  Theuer  is  Commander  and  Director  of  CERL,  and  Dr.  L.  R. 
Shaffer  is  Technical  Director. 
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The  mechanical  properties  (elastic,  plastic,  creep  and  fracture) 
of  cellular  solids  or  foams  are  related  to  the  properties  of  the  cell 
wall  material  and  to  the  cell  geometry.  The  properties  are  well 
described  by  simple  formulae.  Such  materials  occur  widely  in  nature 
and  have  many  potential  engineering  applications. 

1 .  THE  STRUCTURE  OF  CELLULAR  SOLIDS 

Making  foams  is  not  difficult.  Most  polymers  can  be  foamed 
easily,  and  techniques  exist  for  doing  the  same  thing  with  ceramics 
and  glasses.  Even  metals  can  be  formed  into  foams. 

Some  foams  are  almost  isotropic,  meaning  that  their  structure 
and  their  properties  have  no  directionality.  Others  are  anisotropic : 
their  structure  is  axisyrametric  (like  cork)  or  orthotropic  (like 
wood);  and  their  properties  reflect  this.  Man-made  foams  tend  to  be 
almost  isotropic.  Natural  cellular  solids  are  rarely  so;  a  single 
piece  of  cork  or  bone  differs  in  strength  and  stiffness  by  a  factor 
of  2  or  more  along  two  directions  at  right  angles. 

There  is  a  second  distinction.  Some  foams  have  closed  cells 
(like  a  soap  foam):  the  solid  material  is  distributed  in  little 
plates  which  form  the  faces  of  the  cells.  Others  have  open  cells 
(li.*e  a  sponge):  the  solid  material  is  distributed  in  little  columns 
or  beams  which  form  the  cell  edges.  The  mechanical  properties  reflect, 
to  some  extent,  this  distribution.  In  reality,  most  man-made  foams, 
even  those  with  closed  cell  faces,  behave  like  open-celled  foams 
because  surface  tension  draws  much  of  the  solid  material  into  the 
cell  edges  during  manufacture. 

The  most  important  aspect  of  the  structure  is  the  relative 
density,  p/p  ,  where  p  is  the  density  of  the  foam  and  p  that  of 

9  9 

the  solid  of  which  the  foam  is  made.  The  mechanical  properties  of 
foams  depend,  above  all  else,  on  the  relative  density.  It  can  vary 
from  I  to  as  little  as  0.0!  .  The  familiar  foamed  plastics  used 
for  packaging  have  a  relative  density  of  around  0.05  . 
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DEFORMATION  MECHANISMS  AND  IDEALISED  FOAM  STRUCTURES 


(1-4) 

Detailed  studies  of  model  foams  have  identified  four 

deformation  modes:  linear  elasticity,  non-linear  elasticity,  plastic 

collapse  and  various  sorts  of  fracture.  The  studies  used  2-dimensional 

hexagonal  cells,  like  those  in  the  diagrams  shown  later  in  this  paper. 

A  hexagonal  networx  is  a  good  starting  point  because  the  modes  of 

deformation  correspond  to  those  of  3-dimensional  foams,  yet  the 

geometry  is  simple  enough  that  a  complete  analysis  is  practical.  The 

understanding  derived  from  the  2-dimensional  models  has  provided  the 

foundation  for  the  analysis  of  the  more  complicated  problem  presented 

by  real,  3-dimensional  foams^  given  below. 

In  the  analysis,  the  3-dimensional  structure  of  the  foam  must  be 

included.  It  is  idealised,  without  loss  of  physically  important 

features,  in  the  models  shown  in  Fig.l.  The  open-cell  foam  is  modelled 

as  a  cubic  array  of  members  of  length  l  and  square  section  of  side  t 

Adjoining  cells  are  staggered  so  that  members  meet  at  mid-points.  The 

relative  density  of  the  cell,  p/p  (where  p  is  the  density  of  the 

s  s 

cell  wall  material)  is  then  given  by: 


and  the  second  moment  of  the  section  of  a  member  (which  we  need  later) 
is  given  by: 


I 


(2) 


A  closed-cell  foam  is  modelled  similarly.  The  square  struts  are 
replaced  by  square  plates  of  side  l  and  thickness  t  .  Adjoining 
cells  are  again  staggered.  Then: 


1 13 
12 


(3) 


As  pointed  out  earlier,  most  foams  behave  more  like  the  first  model, 
because  surface  tension  concentrates  material  into  the  cell  edges 
during  their  manufacture.  We  will  use  it  as  the  basis  of  the 
calculations  given  below,  which  treat  isotropic  foams.  Refinements, 
results  for  the  second  model  and  the  generalisations  to  non-isotropic 
foams  are  given  elsewhere .  Symbols  are  defined  in  Table  1. 


3. 


MECHANICAL  PROPERTIES 


When  a  foam  is  compressed,  che  stress-strain  curve  shows  3  regions 
(Fig. 2).  At  low  strains,  the  foam  deforms  in  a  linear-elastic  way; 
there  is  then  a  plateau  of  deformation  at  almost  constant  stress;  and 
finally  there  is  a  region  of  densification  as  the  cell  walls  crush  to¬ 
gether.  The  extent  of  each  region  depends  on  relative  density  p/ps  • 
Elastic  foams,  plastic  foams,  and  even  brittle  foams  all  have  3-part 
stress-strain  curves  like  this,  though  the  mechanism  causing  the  plateau 
is  different  in  each  case. 

The  deformation  of  a  regular,  2-dimensional  foam  like  that  sketched 
in  Fig. 3,  can  be  analysed  with  precision.  It  is  far  more  difficult  to 
do  the  same  analysis  for  a  3-dimensional  foam,  because  the  response  is 
an  average  of  that  of  cell  walls  of  random  orientation  in  space,  and 
with  a  distribution  of  length  l  and  section  t2  .  It  is  better  to  use 
the  dimensional  argument  given  below,  and  then  rely  on  experiment  to 
determine  a  single,  unknown,  constant  of  proportionality.  This  method, 
which  we  have  used  to  analyse  many  properties  of  foams,  will  become 
clearer  as  it  is  applied  in  the  following  sections. 


3.1  Linear  Elastic  Properties 


When  a  foam  is  loaded,  the  cells  walls  at  first  bend^2’^’2 
Fig. 3  shows  this  bending  for  the  2-dimensional  model;  the  same  bending 
deformation  occurs  in  3-dimensional  foams.  A  force  F  ,  applied  as 
shown,  causes  the  non-vertical  beams  to  deflect  by  6  ,  which  is 
calculated  from  simple  beam  theory  as: 


6 


C,  F  *3 
12  E  I 


(4) 


s 

Here  Cj  is  a  resolution  factor  which  depends  only  on  the  cell  geometry 
and  Es  is  the  Young's  modulus  of  the  solid  cell  wall  material.  For 
open  cell  foam  of  Fig. I (a)  the  stress  is  proportional  to  F/l 2  ,  the 
strain  to  <5/i  .  The  second  moment  of  area,  I  ,  is  proportional  to  t4 
(eqn.  2)  giving: 


Using  eqn.(l)  for  the  density,  we  find: 


E 

E 

C2<f>2 

s 

s 

(5) 


where  C2  is  a  constant.  The  shear  modulus  scales  in  a  similar  way, 


*  •* .  r.  v-. 


_ 
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because  shear  deformation  in  a  foam  also  causes  simple  bending  of  the 
cell  walls  ^  . 

Data  are  compared  with  eqn.(5)  in  Fig. 4.  It  shows  Young's  moduli 

for  polymeric  and  ceramic  foams  plotted  against  o/pg  on  logarithmic 

scales.  The  full  line  is  a  plot  of  eqn. (5)  with  *  1  ;  it  gives  a 

good  description  of  a  wide  range  of  materials  and  densities.  Note  that 

the  modulus  falls  rapidly  as  the  volume  of  void  space  in  the  foam 

increases.  The  range  of  moduli  practically  available  by  foaming  is 

4 

large:  tt  spans  a  factor  of  10  .  The  modulus  is  important  in  the 

design  of  load-bearing  structures  which  incorporate  foams  (sandwich 
panels  for  example) . 


3.2  Non-Linear  Elastic  Behaviour 


Linear  elasticity,  of  course,  is  limited  to  small  strains, 
typically  5Z  in  compression,  rather  more  in  tension.  Elastomeric 
foams  can  be  compressed  far  beyond  this  point.  The  deformation  is  still 
recoverable  (and  thus  elastic)  but  is  non-linear;  it  is  caused  by  the 
elastic  buckling  of  the  columns  or  plates  which  make  up  the  cell  edges 
or  walls^’^’^’ as  shown  in  Fig. 5,  giving  the  plateau  of  the  stress- 
strain  curve  (Fig. 2).  It  is  exploited  in  cushions  and  packaging  to 
give  a  restoring  force  (which  we  now  calculate)  which  is  independent  of 
displacement. 

The  critical  load  at  which  a  column  of  length  l  ,  Young's  modulus 
Eg  and  second  moment  of  area  I  buckles  is  given  by  Euler's  formula: 
n2  it2  E  I 


(6) 


The  constant  n2  describes  the  degree  of  constraint  at  the  ends  of  the 
column.  If  this  load  is  reached  for  a  layer  of  cells  spanning  the 
section,  they  will  buckle,  initiating  the  elastic  collapse  of  the  foam. 
For  the  3-dimensional  open-cell  foam  of  Fig. I,  the  stress  at 

which  this  occurs  is  proportional  to  F^/i.2  »  so  that: 


E  I 

a*  »  n2  t2  — 
el  *4 


Using  eqns.(l)  and  (2)  we  obtain: 


Several  authors  ’  ’  ’  ’  studying  the  collapse  of  rigid  foams 

sought  to  interpret  their  results  as  elastic  buckling.  Their 
calculations  are  in  the  spirit  of  that  given  here,  but  are 
inappropriate  to  the  materials  they  studied. 
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0*. 

el 

■  m 

E 

C3  ‘f>2 

s 

s 

valid  for  relative  densities  below  0.3  because,  at  higher  densities, 
the  cell  walls  are  too  short  and  stocky  to  buckle. 

Data  for  for  elastomeric  foams  are  compared  with  eqn.(7) 

in  Fig. 6.  They  are  well  fitted  by  the  equation  with  -  0.05  . 

Like  the  modulus,  the  elastic  collapse  stress  spans  a  wide  range:  for 
a  given  material,  a  range  of  10^  is  accessible.  This  is  important 
for  che  design  of  cushions,  padding  and  packaging. 


3. 3  Plastic  Yielding 

Cellular  materials  can  collapse  by  other  mechanisms.  If  the  cell- 
wall  material  is  plastic  (as  are  metals  and  many  polymers)  then  the  foam 
as  a  whole  shows  plastic  behaviour.  It  is  exploited  in  crash  barriers 
and  energy  absorbing  systems. 

Plastic  collapse  occurs  when  the  moment  exerted  on  the  cell  walls 

by  the  force  F  exceeds  the  fully  plastic  moment,  creating  plastic 
.  (2  10  21  22) 

hinges  ’  *  ’  as  shown  in  Fig. 7.  For  a  beam  of  square  section  of 

side  t  ,  che  fully  plastic  moment  is: 


M 


t3 


(8) 


If  a  force  F  acts  with  a  component  normal  to  a  beam  of  length  l  ,  the 
maximum  bending  moment  is  proportional  to  F&  .  The  stress  on  the  foam, 
as  before,  is  proportional  to  F/2.2  .  Combining  these  results  we  find 
the  plastic  collapse  stress  of  the  foam  to  be: 


Using  eqn.(l)  we  obtain: 


C9) 


Data  for  the  plastic  deformation  and  for  the  plastic  indentation 

of  foams  (discussed  next)  are  plotted  in  Fig. 8.  They  are  well  fitted 

by  eqn.(9)  with  C^  -  0.3  for  relative  densities  of  less  than  j  (at 

higher  densities  the  cell  edges  are  too  short  and  stocky  to  bend 

plastically).  The  plastic  collapse  stress  can  be  "tailored",  by  choosing 

3 

a  foam  of  the  right  density,  over  a  range  of  10  It  is  of  primary 
importance  in  the  design  of  crash  padding  and  energy- absorbing  foams. 
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3. 4  Plaatio  Indentation 


Unlike  dense  solids,  which  are  incompressible  when  deformed 

plastically  to  large  strains,  foams  change  their  volume  when 

compressed.  The  cells  of  the  foam  collapse  as  the  foam  is  squeezed, 

so  that  axial  compression  produces  almost  no  lateral  spreading. 

(Poisson's  ratio  for  the  plastic  compression  of  low  density  foam  is, 

typically,  0.04  ( * » 10,24,25) ^  Such  foams  yield  plastically  under  a 

multiaxial  state  of  stress  when  the  maximum  principal  stress,  not  the 

octahedral  shear  stress,  reaches  the  critical  value  a*  calculated  in 

pl 

the  last  section. 

Because  of  this,  the  indentation  hardness  of  a  foam  is  lower  than 
that  of  a  dense  solid  of  the  same  yield  stress:  elements  of  the  foam, 
compressed  beneath  the  indenter,  do  not  expand,  and  so  are  not  constrained 
by  the  surrounding  material  in  the  same  way  that  elements  in  a  dense 
solid  are.  An  analysis  of  the  problem^ shows  that,  for  relative 
densities  less  than  about  0.3  ,  the  intentation  pressure,  or  "hardness" 

H  of  the  foam  is  simply: 


H  -  <Jj1  (10) 

(instead  of  the  result  H  *  3o^  for  a  dense  solid).  Two  experimental 
studies^1**’^  confirm  this  result.  The  data  from  the  second  study  are 
included  in  Fig. 8. 


The  Cruahir 


Brittle  foams  (ceramics,  and  certain  rigid  polymers)  collapse  by 

(26) 

yet  other  mechanisms:  brittle  crushing  in  compression^  ,  brittle 

(27  28) 

fracture  in  tension  ’  .  The  low  crushing  strength  of  refractory 

brick  (a  cellular  solid)  limits  the  loads  than  can  be  applied  to  it; 

and  the  low  fracture  toughness  of  foar  can  cause  problems  when  they 

carry  tensile  loads,  as  they  do  in  w  h  panels. 

* 

Let  the  modulus  of  rupture  ot  >.  -wall  material  be  . 

A  cell  wall  will  then  fail  (Fig. 9)  wh^n  oment  acting  on  it  exceeds: 


(11) 


The  modulus  of  rupture  of  an  elastic  beam,  loaded  in  bending,  is  the 
maximum  surface  stress  in  the  beam  at  the  instant  of  failure.  The 
maximum  stress  is  related  to  the  moment  by  eqn.(il).  The  modulus  of 
rupture  for  a  brittle  solid  is  often  close  to  the  tensile  fracture 
strength. 


14 


As  before,  e  force  F  ,  acting  with  a  component  normal  to  the  wall  of 
length  l  ,  exerts  a  moment  which  is  proportional  to  Ft  .  The  stress 
on  the  foam,  as  before,  is  proportional  to  F/l2  .  Combining  these 
results  we  find  that  collapse  by  crushing  will  occur  at  the  stress: 


3 


Using  eqn. (I),  we  obtain: 


(12) 


There  are  few  experimental  measurements  of  the  crushing  of  brittle 

foams.  The  limited  data^'^,  shown  in  Fig.  10,  are  insufficient  to 

give  confidence  that  eqn. (12)  is  a  good  description.  But  certain  ocher 

observations  do  suggest  that  the  model  has  the  correct  physical 

ingredients.  First  (and  remarkably)  the  tensile  and  compressive 

fracture  stress  of  unnotched  foams  should,  according  to  the  model,  be 

equal  (for  brittle  solids ,  the  crushing  strength  is  roughly  10  times 

(18) 

the  tensile  strength) .  Measurements  on  foamed v  glasses  show  that 
this  is  so.  And,  second,  an  extension  of  the  model  to  describe  crack 
propagation,  awl  to  predict  the  fracture  toughness  of  foams,  describes 
data  well. 


3.6  Ths  Tonsils  Fraature  Stresses 

Compressive  fracture  is  insensitive  to  defects  -  such  things  as 
flaws,  cracks,  or  a  few  exceptionally  large  cells.  But  this  is  not  so 
for  tensile  fracture.  A  completely  unflawed  sample,  it  is  true,  should 
sustain  a  tensile  stress  as  high  as  that  given  by  eqn.  (12)  before  failing. 

But  if  it  contained  a  crack  or  flaw,  then  the  stress  concentration  it  induces 
will  fracture  cell  walls  locally,  extending  the  flaw  and  leading  to  sudden 
fracture  (Fig.  11) . 

When  a  brittle  foam  is  loaded,  the  cell  walls  at  first  deform  elastically. 
The  load  is  transmitted  through  the  foam  as  a  set  of  discrete  forces  and  mom¬ 
ents  acting  on  cell  walls.  But,  since  the  foam  is  linear  elastic  until  the 
cell  walls  buckle,  the  average  force  and  moment  on  a  given  cell  wall  can 
be  calculated  from  the  stress  field  in  the  equivalent  linear -el as tic  con¬ 
tinuum.  We  solve  the  discrete  problem  by  taking  the  solution  of  the 
equivalent  continuum  problem  (just  as  we  do  on  a  smaller  scale  in  replac¬ 
ing  the  discrete  bonds  between  atoms  by  a  continuum)  and  using  it  to 
calculate  the  forces  and  moments  on  the  discrete  cell  walls. 


H>j 


S3 
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A  crack  of  length  a  in  an  elastic  solid,  lying  normal  to  a  remote 

o 

tensile  stress  a  ,  creates  a  singular  stress  field: 


at  a  distance  r  from  its  tip.  Consider  the  first  unbroked  cell  wall, 
which  we  take  to  be  1/2  beyond  the  tip;  it  is  subjected  to  a  force: 


where 


F  «  aZ 2 
o  =  a"  /iZl 


This  exerts  a  bending  moment  on  the  walls  marked  A  and  C,  and  bending 
moment  plus  a  tensile  stress  on  the  wall  marked  B  (Fig.  11).  If,  as 
before,  the  walls  fail  when  the  moment,  proportional  to  F,  exceeds  the 
fracture  moment  given  by  eqn.  (11)  ,  then  the  crack  advances.  Assembling 
these  results  gives: 

»  t3  <*f 

a  *  7T  '  — 

4  /a/i 

where  CTf,  as  before,  is  the  modulus  of  rupture  of  the  cell  wall. 

Using  eqn.  (1)  we  obtain: 

«  o  ^ 

a  s  Cg  (£-)  (14) 

ps  /a / i 

where  a”  is  the  remote  stress  which  will  cause  the  crack  to  propagate. 

The  equation  is  valid  only  when  a  >  £  (otherwise  there  is  no  crack) . 

The  result  is  just  what  would  be  expected.  The  fracture  strength 
equals  that  for  the  unnotched  material  (eqn.  12)  if  the  crack  size  is 
equal  to  the  cell  size.  In  an  open-cell  foam,  a  wall  is  either  broken, 
or  it  is  not.  If  the  crack  size  is  less  than  the  cell  size,  no  cell  walls 
are  broken  and  the  foam  is  undamaged.  Comparing  eqn.  (14)  with  the  definition 
of  Kj£  (the  plane-strain  toughness)  for  a  through-crack  in  an  infinite  sheet 


we  find 


<?£  A r£ 


C8  (S-)' 
K  s 


.  •  « • 


Data  for  crack  propagation  in  a  brittle  polymeric  foam  are  plotted 
in  Fig.  12.  Eqn.  (IS),  with  Cg  =»  0.65,  gives  a  good  description  of 
the  data.  Unlike  all  the  other  properties,  depends  on  cell  size. 

This,  too,  is  expected:  contains  a  dimension  of  length;  the  only 

lengths  which  enter  the  problem  are  the  cell  size  l  and  the  wall 
thickness  t,  but  that  is  related  to  l  through  the  relative  den¬ 
sity  p/pg  (eqn.  1) .  The  result  is  important  whatever  foams  are  used 
in  load-bearing  applications  in  which  tensile  forces  appear  -  as  they 
do  in  sandwich  panels. 


4.  OVERALL  MECHANICAL  RESPONSE :  DEFORMATION  MAPS 

When  an  elastomeric  foam  is  compressed,  it  first  deforms  in  a 
linear-elastic  way;  then  its  cells  buckle  to  give  non-linear  elasticity; 
and,  finally,  the  cells  collpase  completely  and  the  stress  rises 
rapidly  as  opposing  cell  walls  are  forced  together.  A  plastic  foam 
behaves  in  a  somewhat  similar  way,  except  that,  now,  linear  elasticity 
is  followed  by  plastic  collapse,  and,  finally,  the  forcing  together  of 
the  cell  walls.  With  brittle  foams,  progressive  crushing  can  again  lead 
to  a  plateau,  ending  when  the  material  is  completely  crushed.  The 
relevant  formuli  for  open  cell  foams  are  summarised  in  Table  3,  together 
with  the  equivalent  results  for  true  closed-cell  foams  (Fig. lb).  Most 
man-made  foams  behave  mechanically  as  if  they  had  open  cells,  because 
surface  tension  (or  other  factors)  concentrates  the  solid  into  the  cell 
edges . 

The  extent  of  each  phase  of  deformation  depends  on  the  relative 
density.  It  is  convenient  to  display  this  behaviour  as  a  map  with  axes 
compressive  stress  and  compressive  strain,  showing  the  fields  in  which 

each  mechanism  is  dominant.  Superimposed  on  the  field  are  contours  of 
constant  (initial)  relative  density. 

4.1  Deformation  Map  for  Elastomeric  Foams 


Fig. 13  shows  a  map  for  elastomeric  foams.  The  linear  elastic 
regime  terminates  when  elastic  buckling  starts.  The  boundary  of  this 


field  (heavy  solid  line)  lies  at  che  strain  at  which  elastic  collapse 
starts.  Froms  eqns.(5)  and  (7),  this  strain  is: 


e  -  C3  »  0.05  (16) 

At  relative  densities  above  0.3  the  cell  walls  become  so  stocky  that 
they  can  no  longer  buckle  elastically.  The  curvature  of  the  field 
boundary  is  such  as  to  make  the  linear-elastic  loading  line  for 
p/pg  *  0.3  tangent  to  the  boundary. 

The  field  of  elastic  buckling  ends  at  the  strain  at  which  the 
foam  finally  "bottoms  out",  or  densifies,  with  a  rapid  increase  of  load 
with  displacement.  This  starts  when  the  folding  of  the  cells  is  so 
great  that  the  walls  begin  to  touch.  We  find  that  this  begins  when  the 
foam  has  been  compressed  to  a  new  relative  density  of  about  0.5  (that 
is,  the  void  space  occupies  half  the  volume),  and  it  is  complete  when 
the  foam  has  been  compressed  to  a  new  relative  density  of  i  (no  void 
space  left). 

During  elastic  buckling  the  foam  compresses  axially  with  no  lateral 

spreading  (v  ■  0) .  Then  the  relative  density  after  a  nominal 
* 

compressure  strain  e  is  just  (p/pg)  (1/(1 -e) )  .  Equating  this  to 
0.5  gives  the  strain  at  which  densification  starts: 


e 


s 


and  equating  it  to  I  gives  the  strain  at  which  densification  is 
complete: 


(17) 


E  -  1  -  -2.  (18) 

C  ps 

(where  p/p  ,  of  course,  is  the  initial  relative  density).  Within  the 
s 

elastic  buckling  field,  the  stress  is  related  to  the  density  by  eqn.(7); 
using  eqn.(17)  gives  the  equation  of  the  field  boundary  (heavy  solid 
line)  for  the  start  of  densification: 


cs  *  1  "  10 

s 


(19) 


Nominal  compressive  strain  e  ■  (h0-h)/hQ  where  hQ  is  the  original 
height  and  h  the  height  after  a  compressive  strain  of  e  . 


The  contour*  ara  stress-strain  curves  for  foams  of  relative 
-2 

density  ve tween  10  and  I  .  They  show  a  linear  elastic  regime 
(aqn.S),  and  a  plateau  corresponding  to  elastic  buckling  (eqn.7);  they 
start  to  band  upwards  when  densification  starts  (eqn.19)  and  approach 
a  limiting  slope  of  E  when  densification  is  complete  (eqn.18). 

Within  the  field  of  elastic  buckling  the  material  can  exist  in  two 
states  at  almost  the  same  stress  (rather  like  the  p  -  v  response  of 
an  ideal  gas).  It  collapses  in  bands  which  broaden  as  strain  increases: 
Fig. 14,  of  cork,  illustrates  this. 

The  figure  describes  the  overall  response  of  all  isotropic  foams 
in  compression.  Elastomeric  foams  in  tension  show  roughly  linear- 
elastic  response  to  rupture. 


4.2  Plaatia  Foams 

Plastic  foams  (Fig. IS),  like  the  elastic  ones,  show  three  regions: 
linear  elasticity,  plastic  collapse,  and  densification  -  though  now  the 
strain  beyond  the  linear-elastic  regime  is  not  recoverable. 

The  boundary  of  the  linear-elastic  field  (heavy  line)  is  obtained 
from  eqns. (5)  and  (9);  its  equation  is: 

f  -  (0-3/)4-7  (20) 

s  s  e3 

In  constructing  this  map  we  have  taken  ay/Eg  to  be  0.05  .  Next  to 
the  linear-elastic  field  is  the  field  of  plastic  collapse.  As  before, 
two  states  of  strain  coexist  at  almost  the  same  stress,  so  that 
complete  collapse  of  part  of  the  structure  can  occur  while  the  rest  is 
still  elastic  (Fig.  16);  the  bands  of  dense  material  broaden  with 
increasing  strain.  Densification  starts  (as  before)  when  the  cell 
walls  touch  (eqn.17)  and  is  complete  when  the  relative  density  reaches 
1  (eqn.19).  The  field  boundary  (heavy  line)  defining  the  start  of 

densification  is  given,  by  the  arguments  leading  to  eqn.(19),  by: 

c  -  1-4.5 

s  s 

Superimposed  on  the  fields  are  stress-strain  curves  for  foams  of 
initial  relative  density  between  0.05  and  1  .  They  show  a  linear 
elastic  regime  (eqn.5)  and  a  plateau  corresponding  to  plastic  collapse 
(eqn.9);  they  start  to  bend  upward  when  densification  starts  (eqn.21) 
and  reach  a  limiting  slope  of  Eg  when  densification  is  complete 
(eqn. 18) . 


The  figure  shows  Che  overall  response  of  isotropic,  plastic  foams 
in  compression.  It  is  less  general  than  the  map  for  elastomeric  foams 
because  it  must  be  constructed  for  a  particular  value  of  a  /E  .  But 

y  s 

the  equations  show  Chat  the  boundaries  are  not  very  sensitive  to  its 
value,  and,  for  a  given  material,  the  diagram  shows  the  behaviour  for 
all  densities. 

4.3  Rigid  Foams 

Rigid  foams  show  linear-elastic  behaviour  (eqn.5)  to  fracture. 

In  compression,  the  foam  crushes  at  constant  stress  (eqn.12),  and  since 
the  crushing  equation  has  the  same  form  as  that  for  plastic  collapse, 
the  behaviour  will  resemble  that  of  Fig.  17.  If  the  foam  is  contained, 
it  will  densify  at  the  strain  given  approximately  by  eqn.(21),  with 
oy/Eg  replaced  by  cTf/Eg  ■ 

4.4  Selection  of  Foams  for  Mechanical  Design 

The  deformation-mode  maps  combine  the  data  for  a  class  of  cellular 
solid  with  the  understanding  of  the  mechanisms  of  deformation.  Flexible, 
or  elastomeric  foams,  are  widely  used  for  cushions,  padding  and 
packaging.  In  these  applications,  it  is  the  plateau  stress,  (eqn.7) 

and  the  extent  of  the  plateau  (given  by  .inverting  eqn.20)  which  are 
important.  All  this  information  is  summarised,  for  all  flexible  foams, 
in  a  single  map  (Fig. 13).  It  shows  how,  by  judicious  choice  of  E  , 
p  ,  and  p/p  ,  a  foam  best  suited  to  the  application  can  be  chosen. 
Plastic  or  "rigid”  foams  are  used  for  crash  barriers  and  energy  absorbing 
systems:  the  aim  is  to  absorb  energy  while  preventing  the  local  stresses 
or  loads  from  rising  above  the  level  which  will  injure  or  damage  the 
contents  of  the  package.  For  this  purpose,  a  long  plateau  at  a  level 
(eqn.9)  below  the  damage  level,  is  required.  Again  the  necessary 
information  is  sunsnarised  in  the  diagram  (Fig. 15)  which  helps  in  the 
rational  choice  of  foam  material  and  density.  The  crushing  of  brittle 
foams  is  more  of  a  problem  than  a  property  to  be  desired.  The  low 
crushing  strength  of  refractory  bricks  limits  the  load  which  they  can 
carry;  and  the  crushing  of  bone  is  a  problem  when  large  compressive 
loads  are  applied  to  it.  The  deformation-mode  map  for  a  class  of 
brittle  cellular  solids  is  identical  to  that  for  plastic  foams  (Fig. 15) 
with  ay/E3  replaced  by  af/Eg  »  though  the  natural  variability  of 
brittle  solids  means  that  the  predictive  precision  is  lower. 

There  is  considerable  scope  for  developing  the  diagrams  further, 
and  for  investigating  diagrams  to  describe  energy  absorption  by  foams. 
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CONCLUSIONS 


(a)  When  a  cellular  solid  vith  a  relative  density  below  0.30  is 
compressed,  it  shows  a  stress-strain  curve  with  three  parts:  a  linear- 
elastic  part,  a  long  plateau,  and  a  regime  of  final  densification.  The 
properties  of  an  entire  class  of  such  solids  can  be  summarised  as  a 
deformation  mode  map,  which  shows  how  each  of  the  three  parts  changes 
as  the  relative  density  changes.  The  method  can  be  applied  to  natural 
materials,  such  as  woods,  to  give  a  diagram  which  summarises  the 
properties  of  that  class  of  material. 

(b)  The  mechanical  properties  of  cellular  solids  can  be  modelled  with 
precision.  The  models  lead  to  constitutive  laws  (stress-strain 
relations)  which  have  been  thoroughly  tested  in  simple  compression. 

The  response  to  multiaxial  loading  can  be  inferred,  and  the  limited  data 
are  consistent  with  the  inferences. 


(c)  The  constitutive  laws  establish  a  basis  for  mechanical  design  with 
foams.  The  deformation-mode  maps  give  a  compact  summary  of  properties, 
and  give  guidance  in  selecting  foams  for  a  specific  application. 
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TABLE  1:  SYMBOLS  AND  ONITS 

Density  of  foam  (kg/m3) 

Density  of  cell-wall  material  (kg/m3) 

Young 1 s  Modulus  of  foam  (MPa) 

Young's  Modulus  of  cell-wall  material  (MP3) 
Elastic  collapse  stress  of  elastomeric  foam  (MPa) 
Plastic  collapse  stress  of  plastic  foam  (MPa) 
Yield  strength  of  cell-wall  material  (MPa) 

Strain  rata  (s“l) 

Creep  constants  Cs"1,  MPa,  -  ) 

Crushing  strength  of  brittle  foam  (MPa) 

Modulus  of  rupture  of  cell-wall  material  (MPa) 
Fracture  toughness  of  foam  (MPa  m*1) 

Cell  wall  thickness  (m) 

Cell  size  or  cell  wall  length  (m) 

Half-length  of  crack  (a) 

Second  moment  of  area  of  cell  wall  (m4) 

Force  acting  on  a  cell  wall  (H) 

Euler  buckling  load  for  cell  wall  (N) 

Fully  plastic  moment  of  cell  wall  (Nm) 

Moment  which  will  just  fracture  cell  wall  (Nm) 
Dimensionless  constants 
Dimensionless  constants 
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TABLE  3:  THE  EQUATIONS  FOR  STIFFNESS  AND  STRENGTH  OF  FOAMS 


In  most  foams  the  solid  is  concentrated  into  the  cell  edges;  then  the 
open-cell  formulae  are  appropriate.  The  maps  shown  in  Figs. 15  and  17 
were  constructed  using  them. 

+These  equations  refer  to  true,  closed-cell  foams  with  no  thickening  of 
the  cell  edges.  They  are  derived  by  the  method  given  in  the  text,  but 
using  eqn. (3)  in  place  of  eqns.(l)  and  (2).  The  quantities  Bj  to  B 
are  dimensionless  constants. 


Tlx*  3 -dimensional  structure  of  open  and  closed- 
call  foams,  idealised.  The  cell  walls  meet  so  chat 
loads  cause  bending  moments  to  be  applied  to  the 
cell  walls.  Most  foams  behave  like  the  open-cell 
foam. 


The  typical  shape  of  ehe  stress-strain  curve  for  a  foam 


RELATIVE  DENSITY  o/o. 


Th«  relative  Young's  modulus,  S/E  ,  plotted  against 
relative  density,  p/ps<  Open  symiols  represent  open- 
cell  polymer  foams;  shaded  symbols  represent  closed-cell 
polymer  foams,  crosses  are  ceramic  and  glass  foams.  (The  ref¬ 
erences,  in  order,,  are:  1,5,3,11-17;  normalising  constants  in 
Table  2 . ) 


The  relative  elastic  collapse  stress,  <J«l/Es«  plotted 
against  relative  density,  p/ps.  Open  symbols  represent 
open-cell  polymer  foams;  shaded  symbols  represent  closed 

_ 1  e  «  m  m _ 


cell 

norma 


folymer  foams.  (The  references,  in  order,  are:  1,5 
ising  constants  in  Table  2. . ) 


Fig.  6 
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RELATIVE  DENSITY  f/a, 

?ig.  8  Th«  relative  plastic  collapse  stress ,  o£i  /<jv,  plotted 

against  ralativa  density,  o/ps.  Open  symbols  are  for  metal 
foams r  shaded  symbols  are  for  closed-cell  polymer  foams.  (The 

stant^S^able0?*''  “*  l0'3'22'  15,19, 16*,21i  normalising  co 
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The  crushing  of  a  brittle  foam:  the  modulus  of  rupture 
of  the  cell  walls  is  exceeded,  causing  them  to  fracture 


The  relative  crushing  strength,  Of/of,  plotted  against 
relative  density  o/os,  for  brittle  foams.  (The  references 
in  order,  are:  13,22;  normalising  constants  in  Table  1.) 
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Fig.  11.  Mode  1  crack  propagation  in  a  brittle  foam:  the  field 
of  the  crack  subjects  cell  walls  to  a  bending  moment 
and  to  tension.  When  the  breaking  strength  of  the  cell 
wall  is  exceeded,  the  crack  propagates. 
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FRACTURE  TOUGHNESS 
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RELATIVE  DENSITY  p/c. 

Fig.  12.  The  normalised  fracture  toughness  K^/cf  *''*£  plotted 
against c/Ps  ^or  crack  propagation  in  brittle  foams. 

The  references,  in  order,  are:  27,28;  normalising  constants 

in  Table  2.) 
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PART  2 


DEFORMATION  MAPS  FOR  SOLID  POLYMERS 


Summary 


Data  for  Che  viscoelasdc  modulus  can  be  presented  as  diagrams 
which  sumnarise  how  che  modulus  depends  on  temperature  and  time.  The 
diagrams  are  useful  as  a  way  of  classifying  polymers  by  the  mechanical 
response  and  in  summarising  information  about  the  mechanisms  of 
deformation  of  the  polymer,  and  in  presenting  design  data  in  an  easily 
accessible  way. 


1 .  INTRODUCTION 


A  number  of  deformation  regimes  can  be  identified  for  a  given 
class  of  polymers.  A  linear  polymer  such  as  PMMA  shows  4  distinct 
regimes,  each  characterised  by  a  certain  range  of  modulus: 


Ca) 


Glassy  behaviour  at 
of  around  4  GPa  . 


T  <  T 

g 


characterised  by  a  Young's  modulus 


(b) 


Viscoelastic  behaviour  close  to  T^  ,  described  by  a  coupling  of 
the  W  L  F  viscosity  with  a  suitable  elastic  modulus,  leading  to 
a  time-dependent  viscoelastic  modulus  between  I  and  103  MPa  . 


Cc) 


Rubber  elasticity  at  T  >  T  ,  caused  by  entanglements,  following 

g 

Hooke's  law,  with  a  modulus  around  1  MPa  . 


(d)  Viscous  flow  at  temperatures  well  above  the  glass  temperature, 
following  a  simple  viscous  flow  law. 

We  have  made  a  preliminary  examination  of  the  range  of  operation 
of  each  of  these  mechanisms,  and  the  models  which  have  been  developed 
to  describe  them,  with  a  view  to  assembling  the  information  into  maps. 
One  possible  choice  of  axes  are  the  visooelastio  modulus  E  and 
normalised  temperature  T  ,  where: 


p  a 

*  *  "TOT 
e 

TOT 

where  a  is  the  stress,  e  is  the  total  strain  accumulated  in  a 
time  t  and 


?  -  T/T 

g 

where  T  is  the  glass  temperature  (Fig.l).  The  temperature  axis  runs 

-  ®  -3  4 

from  T  *  0  to  T  ■  2  ,  and  the  modulus  axis  from  10  to  10  MPa  . 

This  is  not  the  only  possible  way  of  presenting  the  results,  but  it  is 

one  that  nicely  displays  the  regimes  of  glassy,  viscoelastic  and 

rubbery  and  viscous  behaviour. 


2. 


THE  GLASSY  MODULUS 


Amorphous  polymers  below  their  glass  temperature  have  moduli 
around  4  GN/m2  .  These  low  values  reflect  the  stretching  of  the 
relatively  weak  van  der  Waals  bonds  which  bind  one  chain  to  another,  and 
by  rotation  about  certain  carbon-carbon  bonds  in  the  chain  backbone  such 
as  to  elongate  the  molecule  in  the  direction  of  loading  (Bowden,  1968; 
Yannos  and  Luise,  1982).  The  stretching  and  bending  of  the  carbon- 
carbon  bonds  requires  much  higher  forces,  and  do  not  contribute  in  an 
important  way  to  the  modulus.  If  the  internal  energy  of  the  polymer, 
per  unit  volume,  is  V(e)  ,  then  Young's  modulus  can  be  calculated  from 
the  change  in  V  with  tensile  strain,  e  : 


(1) 


The  intermolecular  van  der  Waals  bonds  arise  from  dipole-dipole 
induction  and  dispersion  forces  and  bind  one  polymer  chain  to  its 
imnediate  neighbour.  It  is  a  short-range  interaction,  well  modelled  by 
a  Lennard  Jones  potential: 

V(e)  =•  V  (i)12  -  2(^)6}  (2) 


where  e  * 


Substituting  and  differentiating,  we  obtain: 


dV  i  1 3  r  7 

H  '  12  vo  -  <tW>  > 


d2V  1  1*+  t  8 

—  -  —12  V  { 1 3  -  7  (77—)  }  (3) 

.7  o  l +e  I+e 

de^  y 

=■  -36  V 

o 

The  second  major  contribution  to  the  elastic  deformation  arises 
from  bond  rotation.  Polarised  infrared  spectroscopy  on  samples  of 
glassy  polycarbonate  (PC),  annealed  below  the  glass  transition 
temperature,  gives  direct  evidence  for  the  rotation  of  segments  of  a 
chain  about  a  carbon-carbon  bond.  The  bond  rotations  are  modelled  by 
a  3-link  deforming  unit  (Bowden,  1968;  Yannos,  1974)  or  "strophon". 
Rotation  of  the  middle  link  causes  an  overall  change  in  length.  For 
rotation  to  take  place,  the  hindrance  caused  by  the  sidegroups  must  be 
overcome.  This  has  been  modelled  by  a  two-fold  symmetric  potential 
(Flory ,  1969): 

V 

V($)  *  {I  -  cos  n  (<p-<p  )}  (4) 

4  o 
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Using  this  idea,  Yannos  and  Luise  (1982)  have  calculated  the  contributions 
of  inter-  and  intra-mo lecular  forces  to  the  glassy  modulus  at  OK  for 
six  different  amorphous  polymers.  The  relative  contribution  of  the  two 
deformation  modes  depends  on  the  length  of  the  strophon  links:  the 
longer  the  Link  the  more  difficult  it  is  for  the  bond  to  rotate,  causing 
intermo lecular  forces  to  become  more  dominant.  This  may  explain  the 
contention  between  Bowden  (1968),  who  proposes  intramolecular  forces  are 
dominant  in  the  modulus  of  a  glassy  polymer  using  polyethylene  (PE)  data 
and  Yannos  (1974)  who  proposes  intermolecular  forces  are  dominant,  based 
on  PC  data.  PC  has  longer  strophon  links  than  PE. 

Although  such  models  allow  an  estimate  of  the  modulus  at  OK, 
they  are  insufficiently  precise  for  our  purposes.  Instead,  we  use 
values  of  Eq  obtained  by  extrapolating  measurements  of  E  at  higher 
temperatures,  to  OK. 

Increasing  the  temperature  has  two  distinct  effects.  First, 
thermal  expansion  increases  the  molecule  separation  and  lowers  the  van 
der  Waals  restoring  forces  and  the  forces  opposing  rotation  about  carbon- 
carbon  bonds:  this  gives  a  slow  drop  in  modulus,  but  does  not  introduce 
a  rate-dependence.  Second,  the  thermal  energy  of  the  molecules  now 
makes  local,  thermally  activated,  rearrangements  possible,  giving  extra 
(time-dependent)  strain,  and  an  associated  drop  in  modulus.  The  first 
effect  (thermal  expansion)  leads  to  a  roughly  linear  decrease  in  modulus 
with  temperature.  Yannos  and  Luise  (1982)  develop  this  idea,  taking 
account  of  thermal  expansion:  increase  in  temperature  increases  the 
interchain  distance,  reducing  the  force  required  for  intermolecular 
deformation.  The  results,  even  for  this  one  contribution,  are  complex, 
but  simplify  (to  within  ±!0Z)  to: 

E  -  Zq  (1  -  o  £■)  (5) 

g 

where 

T 

Ig  dE 

(X  ■  ■»  -  -- 

E  dT 
o 

It  is  reasonable  to  assume  that  the  rotational  contribution  behaves  in 
a  similar  way,  so  that  eqn.(5)  describes  the  effect  of  thermal  expansion 
on  modulus.  Lee  et  al.  (1975)  review  data  for  moduli  of  glassy 
polymers  and  conclude  that  they  are  well  described  by  a  linear 
temperature  dependence  (eqn.5).  Van  Krevelen  (1976)  finds  a  better 
fit  to  the  empirical  relation: 


E 


E  / (1  +  2  T/T  ) 


(6) 


which  reduces  to  eqn.(5)  (with  a  »  2)  for  low  temperatures,  but  which 
gives  a  steeper  drop  near  T 

S 

The  second  contribution  derives  from  the  thermally  activated 
motion  of  side  groups,  or  of  short  segments  of  chains,  past  each  other, 
giving  a  little  extra  strain.  This  thermally  activated  relaxation  is 
difficult  to  model  because,  at  low  temperatures,  only  a  few  sites  with 
the  lowest  activation  energy,  Q  ,  can  move;  at  a  higher  temperature, 
more  sites  with  higher  Q  become  mobile.  This  spectrum  of  activation 
energies  leads  to  cumbersome  equations.  We  adopt  the  view  that 
polymers  below  T  exhibit  secondary  transitions  in  addition  to  the 
main  glass-rubber  transition.  Below  the  glass  transition,  there  remains 
some  limited  freedom  of  rotation  or  sliding  either  of  short  segments  or 
side  groups.  The  energy  involved  in  these  transitions  is  lower  than 
that  for  full  movement,  so  they  occur  at  lower  temperatures.  If  we 
assume,  for  one  such  process,  a  single  activation  energy  Q  ,  then  if 
the  relaxation  time  at  a  standard  temperature  T*  is  x*  ,  then  at  any 
other  temperature  T,t  is  given  by: 


t*  exp  |  (I  -  ±) 


(7) 


The  total  strain  resulting  from-pure  elastic  deformation  plus  a  visco¬ 
elastic  contribution  that  is  due  to  the  transition,  at  constant  stress, 
is: 


o  AE  ,  t,. 

e  (  IT  (  -  exp  T)} 


(8) 


where  E  is  given  by  eqn.(l).  Thus,  the  relaxation  can  be  included, 
provided  four  pieces  of  information  are  known:  the  activation  energy 
Q  ,  the  amplitude  of  the  relaxation  ie/c  and  the  relaxation  time  x* 
at  temperature  T*  .  Further  relaxations  can  be  included  by  adding 
further  terms  like  that  given  above. 


3.  THE  GLASS-RUBBER  (OR  VISCOUS)  TRANSITION 

When  sufficient  thermal  energy  is  available,  the  van  der  Waals 
bonds  holding  the  molecular  chains  of  a  polymer  together  start  to  melt. 
This  enables  segments  of  the  once  elastically-bonded  chains  to  undergo 
large  sliding  movements  relative  to  each  other,  introducing  an  element 
of  viscous  flow  into  the  deformation  of  the  polymer.  This  appears  as 
a  reduced  modulus. 

In  an  amorphous  polymer,  a  single  empirical  function  can  describe 


Che  temperature  dependence  of  the  mechanical  relaxation  near  the  glass 
transition.  Williams,  Landel  and  Ferry  (1955)  describe  the  ratio  of 
the  mechanical  relaxation  time  t  ,  at  temperature  T  to  the  value  t 


at  a  reference  temperature  Tg  (with  suitable  choice  of  Tg  )  by 


-C  (T-T  ) 

.  ,  T  N  1  s 

1o*(t  J  ”  Cl  +  (T-T  ) 
si  s 


(9) 


where  Tg  =  + 50  .  If  is  chosen  as  the  reference  temperature. 


then  we  find 


m  -! 7 .44 (T-T  ) 

lOK(_I_)  s  log  a  - - - §_ 

g^T  )  log  n(T  5  51.6  +  (T_T  ) 

8  8  g 


(10) 


or 


7. 57 (T-T  ) 


n(T)  »  n(T  )  exp- 

o 


51.6+  (T-T  ) 
8 


(ID 


Note  that  eqn.(li)  implies  chat  both  the  viscosity  and  activation 
energy  for  viscous  flow  become  infinite  at  T  »  T  -  51.6°  . 

s 

The  prediction  is  chat  all  molecular  motion  should  become 

completely  frozen  at  T  <  T  -  51.6  .  What  happens,  of  course,  is  that 

8 

new  mechanisms  of  deformation  take  ovec}-for  example  relaxation  of  the 
carbomethoxy  side  chain  at  about  20°C  and  the  aliphatic  ester  group 
at  T  <  1 50°C  for  PMMA. 

Using  the  above  formulation  data  obtained  (for  example)  at 

different  temperatures  but  at  similar  times  can  be  shifted  to  coincide 

on  the  time  3cale  to  produce  a  'master  curve'  at  a  chosen  temperature. 

One  such  'master  curve'  is  produced  by  McLoughlin  and  Tobolsky  (1952) 

from  stress  relaxation  data  of  PMMA  (Figs.l  and  2). 

Other  theories  have  been  developed  to  account  for  the  temperature 

dependence  of  viscosity  in  the  glass  transition  region.  They  include 

barrier  theories  using  potential  wells  and  statistical  theories 

considering  enthalpy  (see  Arridge  1975  for  a  short  review).  But  the 

WLF  equation  has  been  the  most  successful  over  the  temperature  range 

T  to  T  +  100  K  .  Although  it  was  put  forward  as  an  empirical 
8  8 

equation,  it  has  some  physical  basis.  Consider  the  'free  volume'  of  a 
polymer  system.  The  free  volume,  ,  can  be  thought  of  as  the 
difference  between  the  total  volume  V  and  that  occupied  hy  the  atoms 
when  densely  packed,  Vq  .  The  total  volume  V  is  assumed  to  be  that 
of  packed  cylinders  with  a  radius  equal  to  the  van  der  Waals  radii 
around  the  chains  (V  )  ,  plus  that  volume  associated  with  vibrational 


VVA 


motion.  The  free  volume  (V-V  )  is  free  in  the  sense  that  it  can 

o 

redistribute,  continuously,  due  to  thermal  motion.  It  provides  the 

extra  volume  required  for  larger  scale  motion.  Experimentally,  the 

occupied  volume  increases  linearly  with  temperature  whereas  the  total 

volume  shows  a  discontinuity  at  T  -  T  (Fig. 3).  Doolittle  (1951, 

S 

1952)  has  shown  that 

g 

n  *  n*  exp  (-p— )  ( l ; 

rV 

where  n*  is  the  viscosity  at  some  standard  state  and  is  the 

fractional  free  volume  (V-V  /V)  .  Assuming  a  linear  dependence  of 
fractional  free  volume  above  with  temperature  we  find 


f  +  <x,(T-T  ) 
g  f  g 


(at  T  *  T  ) 
g 


where  a ^  is  the  volume  expansion  coefficient,  from  which 


"n(^rr) 

g 


-B/f  (T-T  ) 

g  g 

f  /a,  +  (T-T  ) 
g  f  g 


which  is  equivalent  to  the  WLF  equation  with 


B/2.3  f 


V  “f 


where  f  is  the  fractional  free  volume  at  T  . 
g  g 

We  include  this  information  into  the  maps  in  the  following  way. 
Consider  the  rheological  model  shown  in  Fig. 4.  It  describes  the  bond¬ 
stretching  elasticity  by  the  spring  Ej  ,  the  viscous  sliding  of 
molecules  (as  given  by  the  WLF  equation)  by  the  dashpot  n(T)  and  the 
rubbery  modulus  by  the  spring  E ^  .  The  full  constitutive  law  for  this 
standard  linear  model  is 

E„  E.+E. 

2  _  a  .  a  /  1  2^  /t  ,  s 

e  +  —  e  -  ■=-  +  -  ( — = — )  (14) 

n  E]  n  e, 

At  constant  load  (creep)  the  solution  is 


c  -  a{~  +  ~  (I  -  exp  -  —  t)} 

E1  E2  n 

Then  the  apparent  modulus  o/e  ,  for  loading  time  t  at  temperature  T 


This  correctly  goes  to  Ej  for  short  times  and  low  temperatures, 
and  to  * or  long  times  and  high  temperatures. 

Under  constant  strain  (stress  relaxation)  the  result  is  slightly 
different.  We  have  instead 

■♦S'W  •  nr 

The  solution  is 

Ej  E.+E_ 

0  -  Eie  {1  -  e-^  (I  -  exp  --7T-  t)} 

and  the  apparent  modulus  is 

E.  (E.+E  ) 

E£.s.  ■  E,  11  -  Tr-ag (1  -  «» -  -hr-  '»  (,6> 

As  before,  E£<.  -  Ej  when  time  is  short  or  temperature  is  low;  and 
E£s  ■  EjE2/(Ej+E2)  when  time  is  long  or  temperature  high.  But  the 
transition  between  the  two  is  different. 

We  now  identify  Ej (T)  with  the  glassy  modulus,  E2(T)  with  the 
rubbery  modulus,  and  n(T)  with  the  WLF  equation.  The  apparent  modulus 
is  calculated  and  plotted  using  either  eqn.(15)  or  eqn.(16). 


4.  THE  RUBBERY  MODULUS 

Above  T^  ,  amorphous  polymers  may  exhibit  a  rubbery  modulus 
plateau  of  around  1  MN/m2  .  This  value  is  close  to  that  of  crosslinked 
rubbers  and  arises  in  a  similar  way. 

The  molecular  chains  in  the  polymer  form  a  network  by  curling  and 
twisting  around  each  other  forming  mechanical  entanglements  which  behave 
rather  like  the  chemical  branch  points  of  a  crosslinked  polymer.  These 
are  called  entanglements  and  a  polymer  in  this  region  of  behaviour  is 
characterised  by  an  average  molecular  weight  between  entanglements,  . 
As  the  molecular  weight,  ,  increases  and  so  more  entanglements  take 
place,  Me  is  reduced  and  a  higher  rubbery  modulus  observed. 

To  predict  the  modulus,  the  entanglement  network  is  assumed  to 
act  like  a  crosslinked  network.  When  the  polymer  is  strained,  the 
entanglement  points  constrain  the  free  slippage  of  molecules  relative 
to  each  other,  and  the  molecular  network  tends  to  align,  reducing  its 
entropy.  The  resistance  to  this  ordering  results  in  the  modulus  of  the 
rubbery  region. 


Standard  texts  (Treloar  1958,  Ward  1970)  give  a  derivation  of  the 

rubbery  modulus  En  for  small  strain.  It  is  found  that: 

K 

3  o  RT 


where  o  is  the  density,  T  the  absolute  temperature  and  R  the  gas 
constant.  This  equation  has  been  used  to  calculate  the  rubbery  plateau 
on  the  maps. 


5.  NEWTONIAN-VISCOUS  FLOW 

At  high  temperature  (T  >  (.2  T  )  the  van  der  Waals  bonds  have 
completely  melted,  and  even  the  entanglement  points  slip  relatively 
easily.  At  low  shear  rates  the  polymer  flows  in  a  Newtonian  viscous 
manner . 

Apart  from  the  chemical  dependence  of  viscosity,  which  will  not 
change  for  a  specified  polymer  system,  the  viscosity  is  dependent  on 
molecular  weight  and  the  temperature  T  .  The  dependence  of  the 

viscosity  on  is  given  by  (van  Krevelen,  1976): 


log  n 


log  ncr  +  3.4  log(Mw/Mcr)  (if  Mw  >  Mcr  ) 


log  nQ  -  log  ncr  -  log(Mcr/Mw) 


(if  M  <  M  ) 
w  cr 


The  temperature  dependence  of  viscosity  is  given  by: 


nQ  *  n*  exp(Q/RT)  (1? 

The  WLF  equation  can  be  combined  with  eqn.(!8)  to  give,  for  low  shear 


rates, 


log  n 


I7.44(T-T  ) 

log  ncr  +  3.4  iog(Mw/Mcr>  -  5-—  Tcrh 


Data  are  readily  available  to  describe  Newtonian  flow,  although  large 
variations  can  be  found  between  different  investigations. 


S.2  Data  for  PMMA 


Data  from  van  Krevelen  (1976)  given  in  terms 

of  n  (T)/n  (1.2  T  )  are  converted  for  M  =  100,000  using  eqn.(18), 

cr  cr  g  w 

n  (1.2  T  )  *  I . 26 x 1 0 3  Ns/m2  and  M  =  30,000  ,  also  from  van  Krevelen. 
cr  g  cr 


f’. 
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It  is  found  that  the  activation  energies  Q  (in  eqn.(!9))  are  similar, 
157  kJ/mol  for  Diakon  and  136  kJ/mol  for  van  Krevelen's  data. 

The  value  n*  (eqn.(19))  is  dependent  on  the  molecular  weight 
but  can  be  calculated  given  a  value  of  viscosity  at  a  known  temperature 
T  thus : 


n*  *  (log  ncr(T)  +  3.4  log(Mw/Mcr)  }/exp(j^r) 
The  apparent  modulus  is  then  calculated  from 


6.  CONSTRUCTION  OF  THE  MAPS 

The  maps  are  constructed  from  eqns.(15)  and  (16)  together  with 
eqn.(22).  For  each  temperature,  the  modulus  corresponding  to  a  given 
loading  time  t  is  calculated  and  plotted  as  shown  in  Fig. 5. 
Modifications  of  the  equations  give  the  modulus  associated  with  cyclic 
loading.  We  find  that,  to  give  a  good  fit  to  the  data,  the  time 
constant  in  eqns.(.!5)  and  (16): 


n 

E,  ♦  E„ 


requires  modification,  and  we  are  now  working  on  maps  of  doing  this  in 
a  consistent  way,  and  conducting  a  detailed  comparison  with  experimental 


7.  CONCLUSIONS 

Maps  can  be  constructed  which  summarise  the  viscoelastic  modulus 
of  amorphous  polymers  E(T,t)  for  a  wide  range  of  temperature  T  and 
time  t  .  Several  separate  mechanisms  are  involved:  bond  stretching, 
constrained  molecular  sliding,  rubbery  behaviour,  and  true  viscous  flow. 
The  work  thus  far  must  be  thought  of  as  a  feasibility  study,  which 
indicates  that  the  method  has  potential,  and  is  worth  developing  further 


8. 
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PART  3 


FRACTURE  MAPS  FOR  POLYMER-COMPOS LTES 


■Suimna  ry 

The  interaction  between  fracture  mechanisms  of  a  fibrous  composite 
are  complex.  The  effects  of  certain  intrinsic  material  parameters,  like 
the  cohesive  strength  of  the  fibre-matrix  interface  or  ductility  of  the 
matrix  on  the  micromechanisms  of  fracture  and  toughness  are  unclear. 

One  approach  is  to  link  our  understanding  of  the  failure  processes  with 
theoretical  models  of  fracture  to  predict  the  toughness  of  a  composite. 
However,  the  complex  interactions  between  various  intrinsic  material 
parameters  and  a  set  of  equations  would  not  clearly  be  seen.  Instead, 
we  could  construct  maps  based  on  these  models  to  display  information  on 
fracture  and  toughness  in  a  useful  form.  The  fracture  map  would  have 
two  axes  which  are  labelled  using  any  two  of  the  intrinsic  material 
parameters  that  describe  the  fibre  composite  (Table  1).  The  map  would 
be  divided  into  areas  or  fields,  each  one  depicting  a  particular 
failure  mechanism.  The  boundary  between  one  field  and  another  would 
show  a  change  in  the  dominant  mechanism.  Contour  lines  of  predicted 
toughness  could  be  superimposed  onto  the  map  in  a  manner  useful  for 
design  and  material  selection. 


1 .  INTRODUCTION 

Polymeric  materials  containing  long,  strong  fibres,  fracture  by 
a  number  of  mechanisms.  This  section  describes  the  principal  mechanisms, 
and  the  construction  of  maps  which  show  the  important  composite 
properties  affecting  these  failure  processes  and  the  toughness.  They 
are  constructed  from  data  of  the  properties  of  the  fibre,  matrix  and 
interface,  fitted  to  model-based  equations  of  energy  absorption  which 
describe  the  mechanisms.  Where  experimental  data  is  not  available,  for 
example,  the  misfit  strain,  cq  ,  between  fibre  and  matrix,  we  resort  to 
theory  to  predict  cq  .  Alternatively,  where  related  data  is  available, 
for  example,  the  interlaminar  shear  strength,  we  calculate  values  for 
the  strength  and  toughness  of  a  fibre-matrix  adhesive  bond  using 
fracture  mechanics  relationships. 

The  first  part  of  this  section  describes  certain  fracture 

mechanisms  like  fibre-matrix  decohesion  (debonding),  fibre  fracture  and 

fibre  pull-out  and  shows  how  to  predict  two  important  failure  parameters, 
* 

fibre  debond  length,  ,  and  fibre  pull-out  length,  .  The  second 
part  shows  how  these  predictions  of  fibre  lengths  can  be  combined  with 
other  model-based  equations  to  estimate  the  energies  dissipated  for  a 
number  of  fracture  mechanisms.  The  third  part  shows  how  an  estimation 
of  these  fracture  energy  terms  can  be  displayed  in  the  form  of  a  fracture 

Throughout  this  section,  the  word  'fibre'  is  used  to  mean  either  an 
individual  fibre  or  a  bundle  of  fibres. 


map  co  show  Che  dominant  mechanism  of  Coughening  and  che  predicted 
toughness  for  a  particular  set  of  material  properties.  A  case  study  is 
carried  out  which  shows  how  a  map  can  be  used  to  help  identify  the 
origins  of  hygrothermal  aging  phenomena. 

The  selection  of  maps  presented  is  incomplete.  They  can  be 
divided  into  three  families  of  composite:  glass  fibres,  carbon  fibres 
and  Kevlar  fibres  in  epoxy  resin.  However,  the  method  of  analysing  the 
fracture  behaviour  of  fibrous  composites  and  the  construction  of  a  map 
can  be  applied  to  other  composite  systems,  short  fibres  in  thermoplastic 
matrices,  for  instance. 

1.1  Fracture  Mechanisms 

The  fracture  stress  of  a  brittle  solid  like  glass  or  graphite  is 
not  unique.  It  depends,  for  instance,  on  the  distribution  of  flaw 
size  in  the  sample.  Most  of  the  advanced  fibres,  carbon  for  instance, 
show  a  strength  dependence  on  radius  and  length.  The  location  of  a 
fracture  site  along  the  length  of  an  embedded  fibre  is  also  sensitive 
to  the  rate  at  which  the  load  builds  up  in  the  fibre.  This  will  depend 
upon  the  elastic  properties  of  the  fibre  and  matrix,  modulus  and 
Poisson's  ratio,  for  example,  and  shear  stress  at  the  fibre-matrix 
interface.  In  a  monotonic  tensile  test,  the  fibre  is  first  likely  to 
debond.  The  stress  at  which  it  does  so  is  affected  by  the  strength  and 
toughness  in  shear  of  the  fibre-matrix  interface;  also  fibre  modulus 
and  radius.  A  complicated  relationship  therefore  exists  between  che 
properties  of  the  composite  and  the  debond  stress,  ,  and  fracture 
stress,  Oj  ;  and  between  and  two  fibre  lengths,  the  debond 

length  and  pull-out  length. 

Fibre  debonding  and  fibre  pull-out  are  the  underlying  failure 
processes  which  affect  the  fracture  toughness  of  the  composite.  It  is 
convenient  to  describe  the  failure  processes  in  terms  of  che  energy 
absorbed  in  creating  new  interfacial  surface  and  the  fracture  and 
pulling-out  of  broken  fibres  from  their  matrix  sockets.  There  are 
other  energy  terms  that  contribute  to  toughness;  the  stored  elastic 
strain  energy  of  the  fibre,  released  when  it  snaps,  matrix  deformational 

it 

energy,  and  surface  energies  of  che  fibre  and  macrix.  The  mechanisms 


The  energy  associated  with  plastic  flow  in  the  matrix  and  fracture 
surface  of  the  fibre  and  matrix  is  small  compared  to  the  total  work 
to  fracture  the  composite  and  will  be  ignored  here. 


*■. «r. ■r. •' , - . ■ 


of  fibre  debonding,  fibre  breakage  and  fibre  pull-out  are  therefore 

related  in  complicated  ways. 


1 . 2  Fracture  Toughness  Maps 

It  would  be  useful  to  have  some  way  of  summarising,  for  a  given 
composite  system,  information  about  the  dominance  of  each  mechanism 
and  the  important  material  properties  affecting  toughness.  We  have 
done  this  by  constructing  a  map.  It  is  a  diagram  with  axes  of  any  two 
of  the  composite  parameters  (there  are  at  least  15  such  variables,  see 
Table  1).  Fig. 1  shows  fibre  strength  plotted  against  fibre  modulus  for 
glass  fibres  in  epoxy.  The  map  is  divided  into  fields  which  show  the 
regions  of  fibre  strength  and  fibre  modulus  over  which  each  of  the 
failure  mechanisms  is  dominant.  The  letters  PO  =>  pull-out  energy, 

EL  *  elastic  (fibre)  energy,  and  INT  »  interfacial  surface  energy. 

The  field  boundaries  are  the  loci  of  points  at  which  two 
mechanisms  are  contributing  equally  to  the  toughness  of  the  composite. 
Superimposed  on  the  fields  are  contours  of  constant  tzuohness  (kJ/m2): 
these  show  the  predicted  toughness  of  the  composite  by  summing  the 
contributions  of  each  mechanism.  Also  displayed  are  contours  of  constant 
fibre  debond  length  and  fibre  pull-out  length  (mm),  two  failure  parameters 
useful  in  carrying  out  a  failure  analysis  of  a  broken  specimen.  In 
constructing  such  a  diagram,  it  is  assumed  that  there  are  only  two 
variables,  and  all  of  the  other  material  parameters  have  fixed  values. 

Such  a  map  summarises,  in  a  simple  way,  both  the  experimental  and  model- 
based  understanding  of  the  fracture  of  the  composite  in  question. 

2.  DEBOND  LENGTH  AND  PULL-OUT  LENGTH  EQUATIONS 

In  this  section,  we  develop,  with  a  brief  explanation,  the  fibre 
debond  length  and  pull-out  length  equations  used  later  in  the  construction 
of  the  maps.  First,  we  select  for  the  debonding  and  then  the  pull-out 
of  fibres,  an  equation  which  is  based  on  a  physically  sound  microscopic 
model  of  stress  build-up  in  the  fibre.  While  theory  gives  the  form  of 
the  equation,  we  then  resort  to  experimental  data  or  simple  ways  of 
predicting  the  values  of  the  parameters  which  appear  in  it. 

2. 1  Fibre  Debonding  and  Fibre  Pull-Out 

When  a  unidirectional  fibrous  composite  is  loaded  monotonically, 
the  tensile  stress  on  the  fibre  and  the  shear  stress  at  the  fibre-matrix 
interface  both  increase.  It  is  likely,  because  of  the  low  shear 


strength  of  fibre-matrix  interfaces,  that  the  bond  eventually  fails. 

A  cylindrical  crack  at  the  interface  propagates  along  the  fibre  length. 
The  process  is  called  "debonding".  If  mechanical  interaction  between 
a  debonded  fibre  and  matrix  is  still  possible,  by  fractional  forces, 
for  example,  then  load  transfer  to  the  fibre  will  take  place.  (Matrix 
shrinkage  during  fabrication  and  cooling  may  be  the  origin  of  such 
interaction,  which  sets  up  a  radial  compressive  stress  on  the  fibre.) 
This  produces  a  non-uniform  stress  along  the  length  of  the  debonded 
region.  Because  of  the  variability  of  strength,  the  fibre  fractures 
at  its  weakest  point.  The  eventual  propagation  of  a  macroscopic, 
transverse  crack  in  a  composite  therefore  produces  a  fracture  surface 
showing  protruded  fibres  of  variable  length,  because  of  the  variability 
of  fibre  strength.  This  process  is  called  "pull-out"  (Fig. 2). 

2.2 

The  stress  at  which  a  fibre  first  debonds  depends  upon  fibre 

modulus,  E,  ,  fracture  energy  of  the  interface,  G~  (in  shear),  and 
£ 

fibre  radius,  r^  ^  : 

ad  *  (4EfG2c/rf)  (1) 

The  frictional  stress  between  a  debonded  fibre  and  matrix  socket 
is  related  to  the  compressive  radial  stresses  produced  by  the  matrix 

during  cooling  and  shrinkage,  and  Che  elastic  properties  of  the  fibre 

.  .  (2) 

and  matrix.  The  fibre  stress  can  be  described  by  : 

o(x)  •  o  (I  -  e  ®*)  (2) 

P 

where  x  is  distance  from  the  tip  of  an  interfacial  crack,  o  and 

P 

B  are  functions  of  the  elastic  constants  of  the  fibre  and  the  matrix 
which  take  into  account  Poisson  contraction  of  fibre  under  load: 

o.  -  e  E./v.  and  8  -  2uv-E  /E^r,;  (1  +  v  ) 

p  or  r  tmtt  m 

E  and  v  are  the  modulus  and  Poisson  ratio  of  the  matrix,  v,  is 
mm  r 

Poisson  ratio  of  the  fibre  and  u  is  the  coefficient  of  friction  at 
the  interface,  c  is  proportional  to  the  residual  "misfit-strain", 
eQ  ,  between  fibre  and  matrix,  where  =  5/r^  (  6  is  the  difference 
in  the  radius  of  the  matrix  socket  under  stress-free  conditions,  and 
its  enforced  radius  equal  to  r^  ) .  Eqn.(2)  shows  the  rate  of  stress 

build-up  decreases  as  the  load  on  the  fibre  increases,  o  can 
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therefore  be  thought  of  as  the  maximum  fibre  stress  that  can  be  produced 
by  frictional  loading;  i.e.  when  Poisson  contraction  of  the  fibre  equals 
the  residual  strain  in  the  adjacent  matrix. 

For  a  fibre  chat  has  debonded  by  an  amount  x  ,  the  fibre  stress 
can  therefore  be  written  as: 

a(x)  *  a  -  (a  -  a  )  e  (3) 

P  P<i 


Under  monotonic  loading,  the  interfacial  crack  spreads  until  the 
stress  on  the  fibre  reaches  (Fig. 3).  The  fibre  debond  length  is 

therefore  given  by: 


spreads  along  the  entire  length  of  fibre. 


The  process  of  'bundle  debonding'  can  occur  in  preference  to 
'fibre  debonding'.  A  bundle  of  fibres  can  be  thought  of  as  a  'large 
fibre'  having  the  properties  of  the  bundle.  A  bundle  containing 
matrix  between  the  fibres  has  properties  similar  to  those  of  a  uni¬ 
directional  composite  with  Vf  =  0.80  . 

The  interface  toughness  parameter,  G2C  ,  in  eqn.(I),  is  a  linear 
function  of  two  constituent  material  properties: 

G2c  *  +  (a  -  2rf)G2|/a  (5) 

where  Gj  and  G^  are  the  critical  strain  energy  release  rates  for  the 
fracture  of  interface  and  pure  polymer  respectively.  The  spacing 
between  fibre  centres,  a  ,  for  a  square  packing  array  is: 

a  *  (irr2)^/Vf  (6) 

Vf  is  the  volume  fraction  of  fibres  in  the  bundle. 

For  a  typical  epoxy,  G2  *  500  J/m2  .  We  use  the  following  G; 
values  of  50,  2  and  60  J/m2  for  E-glass,  Kevlar  and  high-strength 
carbon  fibres,  respectively.  The  misfit  strain,  ,  between  bundle 
and  matrix  was  claculated  to  be  about  5%  ;  the  bundle  radius,  ,  was 
of  the  order  of  500  pm  •  Table  5  shows  the  predicted  debond  lengths 
based  on  the  above  values  combined  with  eqns.(l),  (4),  (5)  and  (6). 


«*  / 


Under  monotonic  loading,  an  embedded  brittle  fibre  carrying  a 
variable  load  along  its  length  will  either  snap  at  a  large  flaw  with  a 
low  stress  or  at  a  small  flaw  with  a  high  stress.  The  fibre  shown  in 
Fig. 4  will  actually  fracture  at  point  A  ,  away  from  the  surface  of  a 

large  crack,  and  pull-out.  Its  pull-out  length  depends  upon  the  fibre 

stress  distribution  and  the  flaw  size  distribution  along  its  length. 

Under  stress,  a  ,  the  fraction  P(o)  of  a  population  of  fibres 
that  fail,  in  its  simplest  form,  is: 

P(o)  «  1  -  e (”(a/a0)  )  (7) 

However,  the  probability  of  failure  occurring  in  the  itl1  segment  of  a 
debonded  fibre  is  not  simply  P(a^)  .  It  also  depends  upon  the 
probability  that  a  more  highly  stressed  section  has  not  broken  before 
the  flaw  in  the  i  section  has  caused  the  fibre  to  snap.  This  is  given 
by  the  sum  of  I  -  P(a^)  for  all  j  >  i  ,  i.e.  all  sections  more  highly 
stressed  than  section  i  .  The  relative  probability  of  fracture 

occurring  in  section  i  can  be  given  by  an  integral  form  of  the 

cumulative  probability  function: 

rx  rV2 

P(a(x' ) )  (  1  -  P(a(x"))dx")  dx' 


P (cr(x' ) ) 


l  -  P(o(x  ) )dx")  dx’ 


The  pull-out  length  is  given  by: 

*  -  (£,/ 2)  -  x  (9) 

P  d 

and  F(x')  is  the  cumulative  probability  of  x  being  less  than  x'  . 

Consequently,  the  cumulative  probability  distribution  of  pull-out 

length  being  less  than  or  equal  to  l  is: 

P 

1  -  F| (£./ 21  -4)1  (10 

d  p 

For  composites,  where  ,  no  fibre  debonding  occurs.  Thi 

is  found  in  strongly-bonded  systems  like  carbon  fibres  in  epoxy. 
Individual  fibre  pull-out  is  not  possible  but  'bundle  pull-out'  can 
occur  (we  treat  the  bundle  like  a  large  diameter  fibre  where  a ^  for 
the  bundle  is  less  than  for  the  fibre;  and  for  the  bundle  is 

about  80%  of  the  individual  fibre).  The  nature  of  pull-out  is 


therefore  controlled  by  the  debonding  process.  If  separate  fibres  can 
debond,  they  will  pull-out  individually;  glass  and  KevLar  fibres  would 
be  expected  to  do  this.  The  distribution  of  pull-out  can  be  determined 
from  a  knowledge  of  P(o)  (eqn.7)  and  the  fibre  stress  distribution 
(eqn.3).  Evaluation  of  the  integral  (eqn.8)  can  be  carried  out  using 
an  analytical  method  or  by  a  numerical  technique.  As  a  first 
approximation,  the  average  pull-out  length  for  glass  fibres  and  Kevlar 
fibres  in  epoxy  was  found  to  be  related  to  in  a  simple  way: 


(11) 


(12) 


3.  FRACTURE  ENERGY  (TOUGHNESS)  EQUATIONS 

When  a  composite  ply  fractures,  by  the  propagation  of  a  crack  from 
a  pre-existing  notch  or  hole,  for  example,  a  sequence  of  energy  absorbing 
events  occur  in  a  region  surrounding  the  notch-tip.  For  example,  under 
increasing  load,  the  matrix  cracks,  leaving  intact  fibres  between  the 
faces;  the  fibres  debond  and  new  interfacial  surface  is  created;  fibres 
snap  at  weak  points  and  stored  elastic  energy  is  dissipated;  and  most 
broken  fibres  pull-out  of  matrix  sockets  doing  work  against  friction  at 
the  interface  (Fig. 5). 


3. 1  Inter  facial  Surface  Energy  Equation 

The  surface  energy  of  an  interface  is  small,  perhaps  a  few  J/m2. 
However,  the  area  of  debonded  fibres  and  fibre  bundles  can  be  large  and 
a  significant  amount  of  energy  may  be  required  to  produce  new  surface. 
The  contribution  of  the  fibre  (yj)  and  bundle  (y^)  to  this 
toughening  mechanism  is  additive  and  in  proportion  to  the  total  area 
of  debonded  interface: 


(13) 


The  effective  bundle  fracture  energy  is  given  by  the  weighted  average 
of  y |  and  y2  : 
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Yb  "  l(lTrfYi^  +  (a  “  2rf)Y2[/a 
is  the  work  to  fracture  the  matrix  in  shear 
Fibre  Elastic  Energy  Eauation 


Under  monotonic  load,  the  elastic  strain  energy  stored  in  a 
debonded  fibre  and  matrix  increases.  When  the  fibre  snaps,  some  of  this 
energy  is  released  in  the  form  of  heat,  acoustic  energy,  etc.  An 
estimation  of  the  dissipated  energy  can  be  made  from  a  knowledge  of  the 
scored  energy  in  the  fibre  before  and  after  fracture  (Fig. 6).  The 
dissipated  elastic  energy  is  therefore  the  difference  between  the  initial 
and  final  states  of  stress.  The  elastic  energy  stored  in  the  debonded 
length  of  fibre  immediately  before  the  fibre  snaps  is: 


fV2 

U^  *  2  I  irr^  j  <j(x)  j  2/2E^  dx 

*  n 


where 


a(x)  *  a  -  (a  -or.)  e 
P  P  d 

Assuming  that  after  the  fibre  fractures  the  stress  distribution 
is  linear  (Fig. 6), 

<j(x)  -  all  -  x/  (a  -  a  )  | 

m'  p  d 


f(a  -a.) 

(irrjj/Ef)  P  (a (x) 2dx  -  (irr2/Ef )  (a2 (a  -Oj) /3) 
o  p 


Therefore: 
G  , 


U.  -  U. 
l  f 


3.3  Fibre  Pull-Out  Energy  Eauation 


The  stress  to  pull  a  fibre  out  of  its  socket  embedded  at  a  distance 

x  is : 

a(x)  -  a  (1  -  e  Sx) 

P 

and  the  work  to  pull-out  a  single  fibre  over  a  distance  i  is: 


w  =  irr^o(x)  dx 

If  the  pull-out  lengths  have  a  uniform  probability  distribution,  with 


4.  CONSTRUCTION  OF  THE  MAPS 

The  failure  mechanism  which  is  dominant  (meaning  it  contributes 
most  to  the  total  work  to  fracture  the  composite)  depends  on  the 
properties  of  the  fibre,  matrix  and  interface.  It  is  helpful,  for  a 
given  composite  system,  to  show  the  conditions  under  which  each  mechanism 
dominates  and  to  display  both  the  experimental  data  of  toughness  and 
predictions  of  the  model-based  equations.  Such  a  diagram  or  "map" 
summarises,  in  a  way  which  could  not  be  achieved  by  examination  of  a 
series  of  equations,  a  model-based  understanding  of  the  fracture  of 
composites,  together  with  a  display  of  experimental  data.  Ways  of  doing 
this  were  introduced  briefly  in  Section  1 . 


Method  of  Construction 


First,  data  for  the  material's  properties  are  gathered;  fibre 
strength  and  modulus,  fibre  diameter,  matrix  strength,  modulus  and 
toughness,  interfacial  shear  strength  and  toughness,  and  so  on.  Second, 
we  predict  the  fibre  (and  bundle)  debond  length,  ,  and  pull-out 
lengths,  ip  ;  also  the  frictional  stress  distribution  parameters, 
and  0  ,  from  the  material  properties.  Third,  values  of  and  l  , 

together  with  appropriate  values  of  the  material  properties,  are 
inserted  into  the  model-based  equations  to  predict  the  energy  dissipated 
for  each  micromechanism  of  fracture.  Fourth,  a  computer  constructs  a 
map  by  allowing  any  of  the  material  properties  which  affect  a  fracture 


process,  to  be  varied  m  turn  along  the  two  axes  of  the  map,  keeping 

all  other  material  properties  constant.  The  contours  of  constant 

toughness  are  obtained  by  summing  the  energy-based  equations  to  give 

the  total  work  to  break  the  composite.  Fifth,  predictions  of  i  and 

ip  are  displayed  as  contours  of  constant  length.  Sixth,  any 

experimental  data  plots  of  toughness  and  measurements  of  i,  and  i 

a  p 

are  located  on  the  appropriate  map,  if  required. 

Values  of  the  parameters  used  for  E-glass,  Kevlar  49,  high  modulus 
and  high  strength  carbon  fibres  are  listed  in  Table  2.  A  selection 
of  typical  maps  is  shown  in  Figs. 7- 10. 

4. 2  Comparison  between  Maps 

The  predictions  for  typical  composites,  marked  by  a  triangle,  are 
summarised  in  Table  3.  The  models  predict  reasonable  estimates  of 
toughness,  i^  and  i  ,  for  all  fibre  systems  despite  the  wide 
variation  in  fibre  properties.  In  particular,  the  models  predict  the 
toughness  of  carbon  fibre  composites  to  be  primarily  due  to  the  debonding 
of  bundles  of  fibres,  later  pulled-out  after  fracture;  and  the  toughness 
of  Kevlar  and  glass  fibre  composites  to  be  due  to  the  debonding  of 
individual  fibres,  later  pulled-out  after  snapping.  These  predictions 
are  observed  in  practice.  The  predicted  toughness  of  Kevlar  composites 
is  high,  principally  due  to  fibre  pull-out.  In  contrast,  the  fracture 
toughness  of  high  modulus  carbon  fibre  composites  is  low,  derived  from 
the  debonding  of  fibre  bundles.  A  summary  of  the  predictions  is 
displayed  in  Table  3,  where  the  notation  +  indicates  an  increase  in 
toughness  or  fibre  length  with  increasing  material  property. 


4.  3  A  Case  Study 

Glass-epoxy  laminates  |0/90|  were  fractured  in  3-point  bending 
after  exposure  to  a  hot,  humid  environment  for  up  to  100  h  .  The 
toughness,  G£  ,  was  claculated  and  the  fibre  pull-out  lengths  were 
measured  using  an  optical  microscope.  We  observed  a  fall  in  toughness 
from  27  kJ/m2  to  10  kJ/tn2  ,  approximately.  Fig. 11  displays  the 
experimental  data  on  a  facture  map.  We  believe  that  the  changes  in 
toughness  are  principally  due  to  variations  in  fibre  strength  and 
fibre-matrix  misfit  strain,  which  is  read  from  the  map.  This  implies 
that  eq  initially  drops  due  to  moisture  induced  swelling  of  the 
matrix.  Degradation  of  fibre  strength  by  moisture  attack  causes  the 
large  reduction  in  toughness  with  time. 


5. 


CONCLUSIONS 


The  micromechanisms  of  fracture  and  the  toughness  of  a  fibre 
composite  can  be  displayed  in  the  form  of  a  map.  These  maps  show  how 
the  fracture  behaviour  is  affected  by  changes  in  the  properties  of  the 
composite,  give  guidance  in  designing  a  composite  for  optimum  toughness, 
and  provide  insight  into  the  origins  of  environmental  degradation. 

The  application  of  these  maps  to  predict  the  fracture  stress  of 
composite  laminates  containing  holes  and  notches  is  important  and  needs 
to  be  pursued  further. 
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TABLE  2.  TYPICAL  VALUES  OF  MATERIAL  PROPERTIES  USED 


GLASS 

KEVLAR 

CARBON 

HIGH  MOD. 

HI®  STRENGTH 

Bl  (GPa) 

3S 

76 

220 

140 

E2  (GPa) 

S 

S.S 

9.3 

9.8 

G12  (GPa) 

1.9 

2.1 

S.3 

S.S 

u  12 

0.3 

0.34 

0.34 

0.34 

d!  (GPa) 

0.30 

1.34 

0.39 

1.24 

a2  (GPa) 

ss 

29.6 

60 

66 

1 12  (MPa) 

as 

60 

S! 

83 

Fibre  Radius 
r£  (uffl) 

Bundle  Misfit  Strain 

«b 

E-glass 

7 

S 

Kevlar  49 

6 

4.3 

Carbon 

4 

. J 

2.9 

Fibre  Type 

Gj  (J  m'2) 

Carbon  High  Modulus 

24 

( 

High  Strength 

61 

Kevlar  49 

2 

E-glass 

50 

(Figures  in  brackets  axa  not  contoured  on  the  maos] 


FIBRE  TYPE 

rp£  Cm) 

Fpb  (an) 

2db  Coffl) 

G  (kJ  a  "*) 

DOMINANT 

MECHANISM 

GLASS 

0.21 

("*  0.1) 

3.3 

61 

INTERFACE 

KEVLAR 

0.71 

(y  0.3) 

7.2 

240 

PULL-OUT 

HI®  MODULUS  CARBON 

(0) 

0.09 

2.9 

19 

INTERFACE 

HI®  STRENGTH  CARBON 

(0) 

0.22 

7.7 

67 

INTERFACE 

68 
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Pulled-out  single  fibres 


Fig.  5.  Schematic  of  debonded  and  pulled-out  fibres 
and  bundles  from  a  cracked  matrix. 


6.  Stress  profiles  in  an  unbroken  and  broken  fibre. 
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Figure  10.  Cont'd, 
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BUNDLE  RADIUS  (pm) 


